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We study quantum correlations in an isotropic Ising ring under the effects of a transverse 
magnetic field. After characterizing the behavior of two-spin quantum correlations, we 
extend our analysis to global properties of the ring, using a figure of merit for quantum 
correlations that shows enough sensitivity to reveal the drastic changes in the properties 
of the system at criticality. This opens up the possibility to relate statistical properties 
of quantum many-body systems to suitably tailored measures of quantum correlations 
that capture features going far beyond standard quantum entanglement. 
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1. Introduction 

Quantum many-body systems of interacting spin particles embody a splendid sce- 
nario for the study of quantum statistics, the investigation of fundamental questions 
at the genuine multipartite level and the simulation of non-trivial interaction mod- 
els of generally difficult natural accessibility. Moreover, a very successful decade 
of theoretical efforts has crowned quantum spin systems as promising devices for 
the realization of computation and short-haul communication protocols in on-chip 
solid-state quantum devicef3 2 l 3 l 4 l 

Much of the handiness and interest in analyzing chains and rings of interacting 
quantum spins comes from an extensive body of work performed in the area of 
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so-called exactly solvable models, where analytical and numerical techniques for the 
investigation of statistical properties of a vast class of many-body systems have been 
formalized and put in plac^M3H A nice link between such investigations and a gen- 
uine information theoretical viewpoint comes from the study of figures of merit such 
as entanglement in relation to critical properties of interacting many-body systems 
dragged across a quantum phase transition (QPT^Ml"! 11 ! 12 ! 13 ! 14 ! 15 !!! The incr eas- 
ing interest in the understanding of general quantum correlations' 17 ! 18 ! 19 ! 20 ! 21 ! 22 ] and 
their distinctive features that often depart from those characterizing the latter mo- 
tivate a research aiming at establishing a connection between statistical properties 
of quantum many-body systems and general quantum correlations. Some interest- 
ing steps have been performed along these lined 23 l 24 EH, while the topic still deserves 
a more systematic development. 

We thus present a detailed analysis of the quantum correlations in an Ising 
spin ring in transverse magnetic field l 26 l 27 l 28 l and focus our attention on the global 
properties of the system. In fact, we aim at establishing a formal parallel between 
a recent study by some of us on the genuine non-locality content of the ground 
and thermal state of such a systerrPO and the degree of global quantum correla- 
tions set among the spins by the interaction here at hand. We consider different 
quantifiers of general quantum correlations, each capturing the various aspects of 
such a multifaceted problem. We go beyond the study of mere two-spin quantum 
correlations to consider, on the contrary, their multipartite version. To this task, 
we employ a generalization of quantum discord based on a natural extension of 
mutual information to the many-spin scenario^!. We demonstrate the behavior of 
such correlations as a function of the transverse magnetic field when changing the 
number of spins in the ring. Very interestingly, we show that, while tests for mul- 
tipartite non-locality fail to reveal the modifications occurring within the system 
at criticalitjPS global quantum correlations quite efficiently witness the occurrence 
of a QPT. While a detailed study of such a relation is left for further worlPS we 
believe that our analysis sets the ground for a fertile investigation. 

The remainder of the manuscript is organized as follows: in Sec. [2] we briefly 
review the Ising model of spin one-half particles arranged in a ring configuration and 
under the effects of a transverse magnetic field and briefly discuss the technique 
we use for its diagonalization. In Sec. [3] we recall the definition of bipartite and 
multipartite quantum correlations we are going to use in this paper. Sections [4] and 
[5] are devoted to the study of the behavior of two-spins and global correlations, 
respectively, in the ground state of the quantum spin chain. Finally, in Sec. [6] we 
draw our conclusions. 
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2. The Model 



The model under investigation is an isotropic Ising chain in a transverse magnetic 
field. For N coupled spins, the nearest-neighbor Hamiltonian model reads 

N N 

H = -jY,<®K + l + BY,*n (1) 
n— 1 n— 1 

where J (B) is a dimensionless parameters representing the inter-spin coupling 
strength (the global magnetic field) and {<r^, <r^, <r*} is the set of Pauli spin 
operators for spin n=l, .., N. We assume cyclic boundary conditions such that 
<y l N+ i = a\ (i=x,y, z). The model is exactly solvabld?! 26 ! 32 ! and here we briefly 
recall the methodology used to tackle it. In order to diagonalize the Hamiltonian 
we introduce the raising/lowering operator a — (<r* ± i<7^)/2, and then move into 
a new fermionic picture defined by the Jordan- Wigner transformation^ 

4 = (c„) f = {c„, cl} = s nm , (2) 

where 5 nm is the Kronecker delta. These fermionic variables allow the Hamiltonian 
to be written 

N N 1 

k = - J I>« - + + 2S ^ (gt g„ - -). (3) 

n— 1 n— 1 

A transformation to the momentum representation is then performed using a 
Fourier transform while full diagonalization is achieved by final Bogoliubov trans- 
formation introducing new fermionic operators {bk, b\} (k = —N/2, .., N/2 — 1) and 
giving a free-fermion Hamiltonian 

H s = J2e k b{b k -Y^e k , (4) 

k k 



where e k — \J J 2 + B 2 — 2JB cos 4>k and (f>k — ir{2k + l)/N in the subsector with 
an even number of fermions and a slightly different definition for the odd-number 
case^. The ground state of the system is then the state satisfying the eigenvalue 
equation b\ |GSjy) — Vfc, which gives us 

IGS^v) = (g) (cos ^ |00) fc ,_ fc + i sin ^ \U) kt _ k ) (5) 

with tan d k = (— B + J cos <f> k ) / (J sin 4> k ) and |0)^ fc the state with no (one) 

fermion with momentum <p k . The energy of the ground state is An — — ^2 k e k - 

3. Bipartite and Multipartite Quantum Correlations 

In this Section we introduce the main tools used in our investigation. We first review 
two figures of merit for quantum correlations in bipartite systems: the quantum 
discord and the (ameliorated) measurement-induced disturbance. We then move 
to the multipartite scenario and introduce a measure of multipartite non-classical 
correlations proposed by Rulli and Sarandy in Ref. ^1. 
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3.1. Quantum discord, measurement-induced disturbance and its 
improved version 

Quantum discord is based on the idea of quantifying the discrepancy between 
the generalization to the quantum domain of two expressions for the mutual 
informational. Imagine to have a bipartite system described by the density op- 
erator pab with pa {pb) denoting the reduced state of system A (B). The total 
correlations between A and B are quantified by the quantum generalization of mu- 
tual information 

I{pab) = S(p A ) - S(pa\pb) (6) 

where S(pa) = — TrpA log 2 pa is the von Neumann entropy and S(pa\pb) = 
S(pab) — S(pb) is the conditional entropy. Nevertheless, by using a measurement- 
based approach, a second definition of conditional entropy is possible. The appli- 
cation of local projective measurements on a part of the system projects the total 
system in a different state. In particular if the measurement is described by the 
set of projectors {n^}, the conditional density operator (i.e. the state of the total 
system AB conditioned on the measurement outcome labeled by j) is written as 
p A B\j = {lA®H j B )pAB(lA®n 3 B)/Pj: where pj = Tr{(l A ®~tt J B )p A B} is the probabil- 
ity of outcome j and 1 is the identity operator. One thus define the measurement- 
based conditional entropy S(pab\^ j b ) = J2j PjS(fiA\j) with p A \ 3 = Tv[U B pAB]/Pj 
and finds the alternative version of the quantum mutual information 

J(pab)^S(pa)-S(pab\U b ), (7) 

which is often referred to as one-way classical correlation^. The difference between 
quantum mutual information and classical correlations, minimized over the whole 
set of orthogonal projective measurements on B, defines quantum discord as 

V B ^ A {p AB ) = inf [I(pab) - J(pab)}- (8) 
{nil 

It should be noted that the minimization implied in the definition of D B ^ (pab) 
makes its analytical evaluation very difficult. To date, quasi-closed analytic ex- 
pressions are in fact known only for quite a restricted class of two-spin states^. 
The intrinsic asymmetry of Eq. Q can be lifted by considering the symmetrized 
form V — m&x[D A ~ > ' B (pab) , "D b ^ a {pab)], which is null only on so-called classical- 
classical stated, i.e. density matrices that can be written as J2i 3 ;Pij\h with 
{\i)} and {\j}} single-spin orthonormal sets and pij a joint probability distribution 
for indices (i,j), and is thus a faithful indicator of quantum correlations. 

In Ref.^Luo introduced measurement-induced disturbance (MID) as a different 
quantifier of quantum correlations based on the alterations induced on a quantum 
mechanical state by a measurement process. Under a bilocal complete projective 
measurement {H A <g)H B }, a classical state remains invariant, i.e. pab = H(pab) = 
Sife(nf ® H B )pAB(Hf ® n^). On the other hand, any complete local projective 
measurement and in particular one built on the eigenprojectors {11^ j5 U B k } of the 
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reduced states of A and B, fully decoheres a quantum mechanical state, rendering 
it just a classical statistical distribution of probabilities. The idea behind MID 
is thus to quantify non-classicality by evaluating the difference between p A s and 
H-e(pab) = Eifc( n £, l ( 8 )II l,fc)/ 9 ^s(n^ i «)n| ife ). Quantitatively, MID is defined as 

M(pab) = I(pab) - I(U E ( PAB )). (9) 

Evidently, MID is much easier to compute than quantum discord due to the lack of 
any optimization procedure over the set of projective measurements, and represents 
an upper bound to T>. However, such a lack of optimization can lead to inconsisten- 
cies between the two indicators: MID can be non-null and even maximal on states 
exhibiting zero quantum discord. To remove this inconsistency an improved version 
of MID (AMID) has been propose d 35 ! 36 ! that includes the ab initio optimization 
over any possible set of local projectors on part A and B. Therefore the definition 
of AMID is 

A(pab) = inf[I( PAB ) - I{TI(pab))] (10) 

with II = {Hf ® } as before. The quantitative relation between discord and 
AMID, which is faithful by construction, has been explored and experimentally 
demonstrated using a linear optics setup generating a hyperentangled state of four 
photons^. 



3.2. Global quantum discord 

Here we briefly discuss a measure for the global content of non-classical correlations 
in the state of a multipartite system. By noting that the original definition of 
discord^ can be rewritten in terms of relative entropy^!, the following symmetric 
extension of discord can be considered 



V B ^ A ( PAB ) = S(p AB \\U(p AB )) - S{p A \\Il A (p A )) - S( PB \\U B (p B )) 



(11) 



where S(pi\\p2) — Tr[pi log 2 p\ — p\ log 2 P2] is the relative entropy between states 
Pi and p2- The global quantum discord Q~D{p Al An ), which quantifies multipartite 
non-classical correlations in a system built out of the set of parties {^4 n }, is defined 



inf 

{fbl 



S(p Al ... AN \\fl(p Al ... AN )) 



N 



(12) 



where tlj(p Aj ) = II', p \ II', and fl(p Al ... AN ) = J2k^kPM-A N ^k with Il k = 
^iLi^A an d k denoting the index string Eq. (12), where the infimum 



is taken over all possible multi-local projectors Hj, is always non-negative but its 
maximum value depends on the dimension of the total Hilbert space at hand. 
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4. Two-spin Quantum Correlations in the Ground State of an 
Ising Ring 

In this Section we study the behavior of the two-spin quantum correlations. We 
fix the strength coupling parameter J and investigate the quantum correlations 
shared by pairs of spins "extracted" from the ring against the magnetic field B. 
We also study the behavior of quantum correlations when the ring increases in 
size and pairs made out of non-nearest-neighbor spins are considered. We compare 
the indications provided by the various figures of merit for two-qubit non-classical 
correlations discussed in the previous Section and highlight a series of interesting 
features. 

As here we are focusing our attention to two-spin quantum correlations, we can 
take advantage of the fact that the density matrix of any pair of spins can be 
expressed by means of two-point correlation functions as^ 

Pij = (U+ £ X#*i°®*j 6 )/4 (13) 

a,b— 0,x,y ,z 

with x the two-point correlation matrix with entries xtj = (&? ® &j) and of=l. 
In general, getting the expressions of such correlators is a difficult task due to the 
non-local form of the ground state in the Jordan- Wigner representation^. For com- 
pleteness of presentation, we should however mention that in the thermodynamic 
limit of N — > oo, two-point correlations can be expressed in terms of determinants 
of Toeplitz matrices. The symmetries enjoyed by H are such that the only non-zero 
elements of x are 



Xii-j-s 



G—i G—2 
Go G-i 

G.,_9 G.,_3 



■ G_ s 
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• G_r 

1 + A cos 



Vyy 
ll-\-S 



Gi Go 
G2 Gi 



G s G s _i 
G S G- S , 



G-.5+2 

G-s+3 

Gi 



(14) 



where s is an integer representing the number of sites of the ring separating the 
two elements of the pair being considered. By recasting the Hamiltonian model as 



71+1 



G h = 



we have^l 
1 + A cos ( 



d(f> cos(kcf))- 



dcj> sin(fc(/>) 



sin c 



(15) 



with A = J/B and e = yl + A 2 + 2A cos </>. As long as two-spin states are taken into 
account, the analysis of quantum correlations can be done by directly working on 
the two-spin reduced density matrix pij and using the equations above, for arbitrary 
size of the system. In this work, though, we are interested in global properties of the 
system's state, which in principle require arbitrary multipoint correlation functions, 
not achievable through the apparatus described above, while the considerations to 
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be brought forward in this Section will be easily generalized to arbitrarily sized 
systems, we will consider only finite-length systems to illustrate the key points of 
our study and set useful benchmarks for the analysis presented in Sec. [5j We start 
our analysis by comparing the results coming from the use of the three quantifiers 
introduced above. In particular we consider the quantum correlations contained in 
the reduced state of two neighboring spins as measured by quantum discord, MID 
and AMID. In what follows, without affecting the generality of the analysis, we set 
,7=1 in the model in Eq. ([I]) and leave B as a free coefficient. In fact, the relevant 
parameter in the dynamics under scrutiny is the ratio B/J. 

A remark is due in respect to the computation of some of the figures of merit 
addressed here. It turns out that the reduced density matrix of any pair of spins in 
our problem can be written, regardless of the value of the magnetic field and the 
coupling strength, in the X-like form 



(* *\ 



(16) 



with ★ indicating the only non-zero elements of the density matrix. In this case, 
differently from what occurs for general two-qubit states, a semi-closed analytic 
formula for the evaluation of D3~* % (pij) is available^ (notice that T>^ l (pij) can 
be easily calculated with the very same formula by first applying a Swap gate to 
Pi/). While we point the reader to Ref for full details on this, here it is enough 
to mention that we have used such formula for the calculations reported in our 
work, thoroughly checking the corresponding predictions with an exact numerical 
approach. As for MID, the lack of optimization over the local projective bases makes 
its evaluation straightforward and no further comment is needed in this respect. 
Finally, AMID is calculated exploiting once more the X-like form of the reduced 
two-spin states and the formula found for this task in RefP^ (as for discord, we 
have duly checked the consistency of the analytic predictions with those of a fully 
numerical study). 

In Fig. [T] we present the results corresponding to the case of a ring of N = 6 
spins. At B = 0, the reduced state of any two spins a, ft (with a, /? = 1,2, ..,N), 
obtained as p a ^ = Tr {N , } (\GS N ) (GS N \) (with {A'} = {1, 2, .., N}\{a, p}) is an 
equally-weighted mixture of = (|00) + |11))/V2 and \ip + ) = (|10) + \01))/V2. 
While both discord and AMID predict no quantum correlation in such a class of 
states, M. = 1. One can ascribe such a striking inconsistency of results to the fact 
that, for such a class of states, the reduced single-spin states are proportional to 
the identity operator 1. This implies that the eigenprojections are undetermined, 
exposing the evaluation of MID to a rather rough overestimation. As soon as B>0, 
the balance between such two state components is lost and a proper basis of reduced 
eigenprojectors can indeed be found. Nevertheless, MID quite considerably tops 
both T> and A. This demonstrates that MID fails to capture the genuine content 
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V : M 7 A 




Fig. 1. Quantum correlations shared by nearest neighbors in a chain of N = 6 spins plotted 
against B for J = 1. The plot shows the clear inconsistency between MID and more faithful 
measures of non-classical correlations such as T> and A and the upper-bound to discord embodied 
by AMID. 



of quantum correlations even in such a simple two-qubit state and, as such, does 
not embody a faithful figure of merit for our investigation. On the other hand, A 
provides a much more reliable test, being strictly faithful on classical-classical states 
and only slightly overestimating the symmetrized discord (it should be stressed 
that the differences between such measures rise from the fact that they address two 
different questions, from an operational viewpoint). 

As it will be clarified later on, the main point of our analysis is the global 
content of quantum correlations in the ground state of the Ising ring. Although an 
extension of AMID to the multipartite scenario is certainly possible, and is quite 
naturally entailed in the structure of such a figure of merit, the availability of a 
global quantum discord as a plausible multi-spin indicator and the fact that A only 
constitutes an upper bound to quantum discord push us to consider the latter as 
our key tool for the remainder of our analysis. 

We now consider how the quantum correlation content of the reduced state of 
two neighbors depends on the number of spins constituting the whole chain. In Fig. [2] 
we show how discord in the state of nearest neighbors changes when the number of 
spins in the ring grows. Evidently, regardless of the number of sites in the spin ring, 
the region around B ~ J is special, as the maximum of shared quantum correlations 
can be found within it. From a statistical mechanics viewpoint, the valence of such a 
configuration of parameters is quite understandable as it corresponds to the region 
where the correlation length across the ring diverges. Moreover, using an analysis 
based on the first derivative of concurrence and discord with respect to j j 9 ! 10 ! 38 ^ 
it has been predicted that criticality emerges at B = J. As we will show in the 
next Section, the special nature of this region extends to global general quantum 
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Fig. 2. (Color online) Quantum discord shared by nearest neighbors in a chain of N = 3, 4, 5, 6 
and 8 spins (from top to bottom curve) against the magnetic field B. In the inset we show quantum 
discord evaluated at J/B = 1 against the number of spins in the chain between (i) nearest 
neighbors (blue line/circled dot); (ii) next-nearest neighbors (purple line/squared dot); (iii) third 
nearest-neighbors (yellow line/rhombus dot); (iiii) fourth nearest-neighbors (green line/triangular 
dot) 



correlations. 

The inset of Fig. [2] shows T> in increasingly long rings at a set value of the 
magnetic field [B = J = 1) and for spins separated by a growing number of sites 
(from zero up to three). The larger the separation, the smaller the value of discord, 
although the decrease is rather weak, in stark contrast with what happens to two- 
spin entanglement in the very same spin modeP. Moreover, in Fig. [3] we find that 
the position of the peaks revealed above changes with the site-separatiorP^ 



5. Global Quantum Correlations in the Ground State of an Ising 
Ring 

We now tackle the central part of our study, i.e. the analysis of global quantum 
correlations in the spin model at hand. As anticipated, we make use of the global 
quantum discord proposed in Ref. 30 . When approaching the evaluation of QT>, one 
faces the problem encompassed by the calculation of the relative entropy between 
the ground state \GSn) of the ring and its locally-projected version. For large rings, 
this could be computationally demanding, let alone the necessity of a global opti- 
mization over any possible local projective basis. We thus break the calculation in 
a few intermediate steps and algebraic rearrangements that help us in streamlining 
the evaluation of QT>. 

First, we consider each single-spin projector IF^ as resulting from the appli- 
cation of a rotation R n (9j,(fj) on the projectors onto the eigenbasis of <r*, i.e. 
{Pi = It) (tl ,Po = \i) (i\}n- This leaves us with the observation that the post- 
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measurement state of the ring can be rewritten as 



U(\GS N ) (GS N \)= ]T K({9, tp})V k \GS' N ({0, <p})){GS' N ({9, ip})\ V k VJ({9, <p}), 
fe=l 

(17) 

where \GS' N ({9,<p})){GS' N ({9,<p})\ = K({9,ip})\GS N )(GS N \rt({9,v}) is the ro- 
tated ground state of the ring, lt({9, ip}) is the tensor product of all the rotation 
matrices needed to produce the local projectors, Vk is the tensor product of projec- 
tors over the local a z eigenbasis and we have used the notation {9, ip} as a short-cut 
to indicate the whole set of angles that enter the rotations. Clearly, the value of 
label k determines the combination of a z eigenprojectors {po,Pi} to use in order 



to build up the iV-spin projector Vk- In any case, it is obvious from Eq. (17) that 
the application of each Vk on the rotated ground state picks up a single diagonal 
clement A^ of the latter so that 



V k \GS' N ({9,<p}))(GS' N ({9,<p})\V k =XkA = A fc \k) (k\ (18) 

with \k) the TV-spin eigenstate of ®n=i&n determined by our choice of k. 
The set of Afc's, which is straightforward to determine even for a large den- 
sity matrix, embodies the eigenspectrum of II(|G5jv) (GSn\). The latter is ob- 
viously diagonal in the rotated basis {\k({9,ip})) — H({6, <p}) \k)}. By ar- 
ranging such vectors in columns, we form the passage matrix T such that 




Fig. 3. (Color online) Quantum discord as a function of the magnetic field between nearest 
neighbors (red solid line), next-nearest neighbors (blue dashed line), third nearest-neighbors (green 
dotted line) in a chain of 8 spins with J = 1. Inset: Magnification of the main panel in the region 
of B G [0,2]. 
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(b) 



(c) 




2.0 
1.0 



B/J 



3 5 

B/J 



Fig. 4. Global quantum discord against the ratio J/B for Ising rings of increasing size. We have 
taken N = 3 [panel (a)], 4 [panel (b)] and 5 [panel (a)]. The QT> curve at assigned N is a 
universal function of B/J, peaking in a region around the critical point of the transverse model 
we are considering. 

TU(\GS N ) (GS N | )Tt=A=diag[A 0! A x , .., A 2 iv_ 1 ], which allows us to write 



S(\GS N ) (GS N \ \\U(\GS N ) (GS N \))=J2 M*l°g 2 Wfe-T* [(T+ \GS N ) (GS N \ T) log 2 A] 

fc=l 

(19) 

with /ik the k th eigenvalue of \GSn) (GSn\- Despite its innocence, this expression 
simplifies the evaluation of the global quantum discord. Needless to say, we still 
face the necessity for the optimization implied in the definition of QT>. However, the 
evaluation of its expression prior to this step is now reduced to a computationally 
non-demanding problem, and this opens up the possibility to explore the thermal- 
state scenario where the Ising ring is not prepared in its ground state but is affected 
by the influences of a non-zero temperature^. 

In what follows, we restrict our attention to the ground-state case and investi- 
gate the behavior oiQT> against the magnetic term B. The results of our calculations 
for A" = 3, 4, 5 spins are shown in Fig. [3] where we have studied QT> against the 
ratio B/J so as to investigate universality features of such figure of merit with 
respect to this parameter. At small values of B/J, the ground state of the ring 
is locally equivalent to an A^-spin GHZ state^, for which C/D=lPH By increasing 
the ratio B/J, the system undergoes a profound change in the way correlations 
(not just entanglement) are shared, as it goes from such a genuinely multipartite 
entangled state to a fully separable one achieved for B/J — > oo. The global dis- 
cord spectacularly captures the special changes in the correlations among the spins 
occurring close to J=B and signals it with a singularity. The functional form of 
QT> is independent of the size of the ring, although the decay of collective quantum 
correlations is faster for rings of smaller size. Due to the handiness of the analytical 
results described above, we have been able to numerically sample the behavior of 
QT> for A~ up to 8, finding results in qualitative agreement with the analysis above. 
The sharp peak in the proximity of B/J = 1 suggests the existence of a singularity 
of QV at criticality. In turn, this would imply a discontinuous first derivative, thus 
suggesting a second order phase transition, in agreement with what is known for 



January 17, 2012 1:8 WSPC/INSTRUCTION FILE CampbellMazzolaPa- 
ternostroIJQI 



12 Steve Campbell, Laura Mazzola, and Mauro Paternostro 

(a) (b) 




Fig. 5. (Color online) (a) We study the scaling law regulating the maximum value of global 
discord maxC7X> versus the increasing dimension of the spin ring. The dots show the numerical 
results of our calculations, while the line shows the best linear fit compatible with the constraint 
max C/2?=2? ++ =l for N=2. See the body of the paper for details, (b) Value of the magnetic field 
(B/J)* at which max QT> is achieved, plotted against N. 



the Ising model. The trend followed by global discord contrasts quite evidently with 
what has been found in terms of multipartite non-locality^ the degree of violation 
of multipartite Bell-like inequalities is a monotonically decreasing function of the 
magnetic field, thus signaling the profound differences between general quantum 
correlations and non-local ones. 

We have further analyzed the features of QV by looking for some qualitative 
indications on the scaling law regulating the growth of the peak of global discord in 
the region around B/J = 1. By imposing the constraint m&xQT>=T>^=l at N = 2, 
the numerical results of our calculations appear to be well fitted by the linear 
function 

max gV = m(N - 2) + 1 (20) 

with m the slope of the line. We have thus found the best-fit line for the numerical 
data associated with N = 2,.., 7 provided in Table [lj which gives m=0. 693461 [a 
comparison between the best-fit function and the numerical data points is given in 
Fig. [5] (a)]. On the other hand, the finite-size nature of the examples worked out 
here induces some deviations of the value (B/J)* of the magnetic field at which 
the singularity of QT> occurs from the expected criticality point. Such deviations 
reduce with N growing [cf. Fig. [5] (b)] as a power-law behavior. 

In order to build a parallel with the analysis and the similarities between two- 
spin entanglement and discord, it is interesting to compare the behavior of such 

Table 1. Dimension N of the spin ring and associated value of m&xQT) 



TV 2 3 4 5 6 7 



maxgi) 1 1.8296 2.4360 3.0879 3.7095 4.501 
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a global measure of quantum correlations with the predictions of a multipartite 
entanglement measure. To do this, we need a quantum statistical indicator that is 
able to capture the occurrence of a quantum phase transition. To this purpose, we 
adopt the measure for pure states introduced in RefP^. Given a composite system 
that we divide in the bipartition (A,B), we define the measure of entanglement 

S = -log 2 [Trp A ] (21) 

with pa = Txb |^),4b (^| the reduced state of subsystem A {\^) AB is a pure state 
of the whole system). As 1/Trp^ is viewed as the number of terms entering the 
Schmidt decomposition of l^)^, S is operatively interpreted as the number of 
spins effectively entangled (relatively to the considered bipartitionj^H. For a com- 
pletely separable state Trp\= 1, regardless of the bipartition, so that £ =0. The 
general scenario will see such an indicator depend on the choice of (A, B) and it 



is intuitive to consider the statistical average £ of the values achieved by Eq. (21 1 
achieved by exploring all the possible bi-splitting of a system as a measure of the 
multipartite entanglement shared by its elements. Being this measure statistical in 
nature, the predictions arising from it would be strongly dependent also on higher 
statistical moments. The variance of the distribution of entanglement across any 
possible bi-splitting, in this respect, provides an indication of the sensitivity of the 
entanglement-sharing to the particular bipartition being taken. A large variance 
signifies a strong dependence of £ on the choice of splitting. 

When the apparatus described above is applied to the state of the spin-ring sys- 
tem, the typical behavior of £ and the normalized variance AS = A£/(max AS) is 
displayed in Fig. [6] (a) (the specific case shown in the figure is relative to N = 6). For 
vanishing values of B / J, S~l with a zero-width distribution of entanglement across 
the possible bipartitions, a situation that clearly witnesses the GHZ nature of the 
spin-ring state (whose bipartite entanglement is uniformly distributed across the 
various bi-splitting and this achieves AS=0). Equally expectedly, for large values of 
Bj J no entanglement should be found in the system, regardless of the bipartition. 
Again, this is well signalled by our analysis, which gives a vanishing S with a quickly 
decaying associated variance, as the system tends to a fully separable state. The 
intermediate region of values of Bj J is the most interesting one: while S smoothly 
decreases without exhibiting any special behavior (except changing from concave 
to convex), its variance peaks in the region where B~J [cfr. the vertical marker 
located at J = B in Fig. [6] (a)]. Physically, this signifies an increased sensitivity 
of the entanglement-sharing structure to the specific bi-splitting we consider, thus 
marking the quick departure from the regular entanglement distribution typical of 
a GHZ state. Going from an ordered phase (the GHZ one) to a differently ordered 
one (corresponding to fully separable states) as B increases, the system is forced to 
readjust its inner correlations, breaking the symmetry of the way quantum correla- 
tions are distributed among its parties and thus increasing the associated variance. 
This analysis reinforces the idea highlighted above in terms of the behavior of global 
discord, that such a phase readjustment is a genuinely global phenomenon that is 
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(a) (b) (c) 




Fig. 6. (a) Partition-averaged entanglement £ [solid line] and its normalised variance AS [dashed 
one] for a ring of six spins plotted against the ratio B/J. The vertical line marks the point B=J. 
(b) Comparison among the variances AS (as functions of B/J) of the entanglement distributions 
for a ring of size N = 4, 5,6 (solid, dashed and dotted line, respectively). We mark the position 
of the maximum of each curve with a dashed vertical line, (c) Comparison between QT> and AS. 
As B/J increases, we move along the curve as shown by the arrows. 

very well captured by an equally global indicator of quantum correlations. 

We should notice that the amplitude of the curve describing A£ is a function of 
the ring size N, and so is the location of max AS , as it is illustrated in Fig. [6] (b) 
for N — 4, 5, 6. In analogy with the results displayed in Fig. [5] (a), the maxima are 
well fitted by a linear function of N and the shifts in their position decreases with 
N. Such similarities call loud for a more direct comparison between QT> and A£, 
which is shown in Fig.[6](c), where the two figures of merit are plotted against each 
other (for B/J £ [10~ 6 , 6], growing as indicated by the sense of the arrows long the 
curve and, for easiness of calculation, N = 4). The corresponding open hysteresis 
path shows very clearly the sharpness of the peak of QT> and a small (finite-size 
induced) mismatch between the positions of the maxima of global discord and A£. 
Such relative shift decreases as N grows, thus signalling an increasing accuracy in 
determining the actual position of the critical point as the size of the ring increases. 

On one hand, the study above confirms the predictions coming from the use 
of global discord. On the other hand, it is remarkable in showing that differently 
from quantum discord, which is very efficient in pin-pointing the structural phase 
transition of quantum correlations in the system, the actual quantifier of entangle- 
ment chosen here is unable to do so (in analogy with what is found running tests 
of genuinely multipartite non-locality^) . A more refined statistical analysis is nec- 
essary for this task, as shown by the success achieved in using the variance of the 
entanglement distribution. 



6. Conclusions and Further Developments 

We have studied the behavior of general quantum correlations shared by the ele- 
ments of a spin system governed by a transverse Ising model under the influences of 
a collective magnetic field. Our approach was multifaceted: on one hand, we aimed 
at showing that some degree of care is required in the choice of the indicator of non- 
classicality in such a model. Naive choices dictated by the easiness of calculation 
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implied by basis-dependent figures of merit may well lead to misleading physical 
conclusions. On the other hand, we have embarked in a quantitative investigation 
on the content of multipartite quantum correlations across the spins of a given ring 
prepared in its collective ground state, finding that the deep structural changes oc- 
curring close to the ring's critical point are well revealed by a multipartite extension 
of quantum discord. Such conclusions have been confirmed by a statistical analysis 
of the way entanglement is distributed across a ring of a set size: the variance of such 
distribution agrees in an excellent way with the predictions of global discord, while 
the average degree of entanglement is basically oblivious to criticality, similarly to 
the behavior of non-locality indicators for the very same system^!. Our analysis 
paves the way to a more extensive study addressing thermally-affected states and 
aiming at relating the behavior of global discord to the most intimate statistical 
properties of the spin system studied here^. 
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